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We present a universal relation between the flow of a Renyi entropy and the full counting statistics 
of energy transfers. We prove the exact relation for a flow to a system in thermal equilibrium 
that is weakly coupled to an arbitrary time-dependent and non-equilibrium system. The exact 
correspondence, given by this relation, provides a simple protocol to quantify the flows of Shannon 
and Renyi entropies from the measurements of energy transfer statistics. 


Exact correspondences between seemingly different 
concepts play important role in all fields of physics. An 
example is the fluctuation-dissipation theorem, which 
states that the linear response of a system to exter¬ 
nally applied forces corresponds to the system fluctua¬ 
tions m- In the last decade, the fluctuation-dissipation 
theorem has initiated important developments in quan¬ 
tum transport, quantum computation, and other simi¬ 
lar phenomenological theories [3]. This theorem can be 
extended to nonlinear responses [4j and to full counting 
statistics (FCS) [5], giving more extended sets of such re¬ 
lations similar to Crooks’ formula [6]. In this paper we 
present a relation similar to the fluctuation-dissipation 
theorem that provides an exact correspondence between 
the flows of Renyi entropy and FCS of energy transfers. 

In transport theory, stationary flow of a physical quan¬ 
tity can take place from a system into an infinitely large 
system. In the case a quantity is locally conserved in 
each system, its flow is determined only by interaction 
between the two systems [7. The traditional exam¬ 
ples include electric current, which is the flow of charge, 
and energy flow. Moreover, there are other conserved 
quantities which are not physical in a strict sense. An 
example is the generalization of entropy by Renyi into 
Sm = > with p n being the probability to be in 

state n and arbitrary M > 0, [8]. Quantum general¬ 
ization of the Renyi entropy is obviously conserved in a 
system under Hamiltonian evolution, with the Hamilto¬ 
nian involving only the degrees of freedom of this system 
m ■ For a system in thermal equilibrium at temperature 
T this entropy corresponds to the difference of free ener¬ 
gies, i.e., ln^M = F(T) — F(T/M). In non-equilibrium 
thermodynamics, the Renyi entropies have already been 
considered j9]. They have been studied in strongly inter¬ 
acting systems mm 2 ], in particular spin chains mm- 

The Renyi entropies in quantum physics are consid¬ 
ered unphysical, or non-observable, due to their non¬ 
linear dependence on density matrix. So is the Shan¬ 
non entropy, which is derived from the Renyi entropy 
S = liniM^i OSm/OM, 10 . Such quantities cannot be 
determined from immediate measurements; instead their 
quantification seems to be equivalent to determining the 
density matrix. This requires reinitialization of the den¬ 
sity matrix between many successive measurements m • 
Therefore the flows of Renyi entropy between systems 


J~m = — dlnSu/dt are the conserved measures of non¬ 
physical quantities. The same pertains to Shannon en¬ 
tropy flow m- An interesting and non-trivial question is: 
Is there any relation between the flows of Renyi entropy 
and the physical flows? An idea of such relation was first 
put forward by Levitov and Klich in P!> where they 
proposed that the Shannon entropy flow can be quan¬ 
tified from the measurement of full counting statistics 
(FCS) of charge transfers. The validity of this relation 
is restricted to vanishing temperature and obviously to 
the systems where interaction occurs by means of charge 
transfer. This has been further elaborated in El- 

In this paper we present a relation which is similar 
in spirit. It gives a correspondence between the flows 
of Renyi and Shannon entropies and the FCS of energy 
transfer in the limit of weak coupling. From analysis 
of the previous results for entropy production in quan¬ 
tum point contact [T6| [T8 ] and more general perturbative 
derivations in [TO] , the Shannon entropy flow is known to 
be proportional to heat flow in the absence of external 
forces. This is violated in higher-order perturbation se¬ 
ries. This implies that the exact correspondence does not 
hold for strong coupling limit. 


I. DEFINITIONS AND RESULT 

We consider two quantum systems X and Y. We as¬ 
sume that the system X is infinitely large and is kept in 
thermal equilibrium at temperature T. The system Y is 
arbitrary: it can encompass several degrees of freedom 
as well as infinitely many of those. It does not have to 
be in thermal equilibrium and in general is subject to 
time-dependent forces. It is convenient to assume that 
these forces are periodic with period r. However this pe¬ 
riod does not enter explicitly in formulation of our result, 
which is also valid for aperiodic forces. The only require¬ 
ment is that there is a stationary limit of the flows of 
physical quantities to the system X. The stationary limit 
is defined by averaging the instant flow over the period 
r. For aperiodic forces it is determined by averaging over 
sufficiently long time interval. 

The energy transfer is statistical. The FCS of en¬ 
ergy transfers concentrates on the probability P(E tr ,T ) 
to have energy transfer of E tr during time interval T, 
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mm- In the low frequency limit of long T all statis¬ 
tical cumulants of the energy transfer are proportional 
to T and are determined from the generating function 

F(0 = JdE tr P(E tr: T)eM^E tr ) ~ exp(-r/(0). 

The parameter £ is a characteristic parameter and cu¬ 
mulants are given by expansion of /(£) in £ at £ = 0. 

For quantification of the Renyi entropy flow we need 
to define an auxiliary FCS of energy transfer. The most 
general interaction Hamiltonian is H = X n Y n with 
X n being operators in the space of the system in ther¬ 
mal equilibrium, and Y n being those in the space of the 
arbitrary system. Let us replace Y n with average values 
Y n —>• (Y n ). The result in the Hamiltonian is that of 
the equilibrium system subject to time dependent exter¬ 
nal forces. Those induce energy transfers to the system 
to be characterized by a FCS. We discuss below possi¬ 
ble physical realization of the scheme. So we have two 
FCSs. We denote their generating functions with /i(£) 
(incoherent) and / c (£) (coherent). 

Our main result is the following exact correspondence: 

r = im-r « 

which indicates that the Renyi entropy flow of the or¬ 
der M to the system kept at temperature T = l/ksP is 
exactly equal to the difference of FCS of incoherent and 
coherent energy transfers to the system kept at tempera¬ 
ture T/M at the fixed characteristic parameter £*. This 
relation is valid in the limit of weak coupling, where the 
interaction between the systems can be treated pertur- 
batively. 

There is an obvious classical limit of the arbitrary sys¬ 
tem: all operators Y n are just numbers corresponding to 
classical forces acting on the system in thermal equilib¬ 
rium. In this case the dynamics of the system is governed 
by the Hamiltonian in degrees of freedom of the system 
and therefore will be unitary. In this case, the trace of 
any power of density matrix, as is used in the definition 
of Renyi entropy, will not change in time: There will be 
no entropy flow. This result can also be understood from 
the correspondence ([!]): in this case fi = fe¬ 


ll. DERIVATION OF THE RESULT 

Here we discuss the proof of the exact correspondence 
in Eq. 0 in the weak coupling regime where we can 
restrict ourselves to the first non-vanishing order of per¬ 
turbation theory. To start with, we determine the FCS 
generating function of energy transfers using a diagram¬ 
matic representation of a pseudo-density matrix. Then 
we obtain the Renyi entropy flow from a multi-contour 
technique, and we demonstrate the correspondence of the 
two. The general formalism is illustrated by applications 
to two particular types of systems: the simplest quantum 
heat engine and a harmonic oscillator system coupled to 
environments. 


A. Full counting statistics 

Interactions between two systems influences the statis¬ 
tics of conserved quantities such as current and energy 
flows. In our consideration of FCS of energy transfer, 
we follow the lines of reference mana. We specify it to 
our situation where the interaction Hamiltonian between 
system X in thermal equilibrium and an arbitrary sys¬ 
tem Y is given by H = A n T n . The FCS of energy 
transfer in system X during the time interval [0, T\ can 
be determined from the following generating function: 


F r (£)sTY x/ 5 x (T), (2) 

using the dynamics of the pseudo-density matrix p\ 
p x (T) = Tr Y | ^e _i ^» Tdt2 ^ m ( t2_ I)Ymfc)j x 

p( 0) + (3) 

where (3?) denotes (anti-) time order operator. This 
quantity can be rewritten as a Keldysh partition function 
with integral taken over a Keldysh contour. 

Fig. 0 shows that there are four possible diagrams for 
the evolution of p in the second order. Time moves for¬ 
ward from left to right. Each diagram contains a double¬ 
contour, the outer (inner) represents the evolution of sys¬ 
tem X (K). Determining FCS for system X requires to 
shift X operators in time with ±£/2, with £ being the 
characteristic parameter. The value of the shift is oppo¬ 
site for forward and backward contours corresponding to 
time evolution of bra and ket states. 

Let us consider second order perturbation for dp x /dt. 
An element of this diagram is the average (X(t)X(t f )). 
This average is performed over the states of thermal equi¬ 
librium. We define 

(4) 

The spectral density is Smn{u>) = J d(t — 

t') exp(ico(t — t'))Smn (t — t'). Since t — t' is large 
in Eq. 0 we can shift the lower bound 0 — oo. Due 
to Markov approximation we can replace p x (0) with 

pAT). 

This allows us to compute the mean value of generating 
function /(£) = —Tr(dp x (t)/dt)/Trp x (t) averaged over 
the period r. This can be explicitly obtained from Eqs. 
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FIG. 1: Diagrammatic illustration of dynamics of p using 
multi-contour evolution. Arrows indicate the direction of fol¬ 
lowing Keldysh contour. In each diagram two interactions 
take place at vertices m and n and their correlator link the 
vertices by a solid line. Auxiliary times d=£/2 shift the times 
when interactions act on the contour of X. 


([2| and ^ and the diagrams in Fig. ([!]): 


m = ~- [dtf my: 

J 0 J —oo 



x (Y m {t')Y n (t) p Y ) 


+ 




X m (t+l)X n (t' + ^)p 


x(Y m (t)Y n (t 




( 5 ) 


with the definition 


y mn (w) = - [ dt [ dt' 

T Jo J — oo 

| (Y m (f) Y n (tf e-W-V + (t m (t) Y n (0) J 


( 7 ) 


FCS of coherent energy transfer : Since a driving force 
is externally applied, another type of energy exchange 
is possible to take place between the driving force and 
the system X. In this sense, two opposite energy trans¬ 
fers occur between the external force and system X, one 
at t and the other at t' . The two transfers are corre¬ 
lated from within the system X. One way to consider 
this energy transfer is to replace Y n with the driving 
energy: Y n -D (Y n ). The interaction Hamiltonian is 
H = X n (Y n ) +h.c. The full counting statistics of co¬ 
herent energy transfers can be described in a similar way 
as that of incoherent energy transfer discussed above: 

IF (0 = - E / ^ Y iui - !) S%1 M Y m„M, 

ran J 

(8) 



| (Y m {t')) (Y n {t)) e-^-^ 

+ (y m (t))(f„(0)e iu ’ ( ‘- t ' ) j (9) 


with Ymn being spectral density of the forces acting on 
the system X. 

Using the relation between S and response function 
Xmn (see Appendix [B]) one can rewrite the coherent FCS 
in more comprehensive way: 

fF(0 = ~ j ^E^nH Y mnHx 

mn 

- 1) n (, -jj/T ) + (e iw « - l) (n (w/T) + 1)] 

( 10 ) 


Statistical cumulants C n can be determined from the 
FCS generating functions from C n = i n d n //d£ n at f = 0. 


FCS of incoherent energy transfer : The incoherent en¬ 
ergy transfer takes place between the system X and Y 
using the interaction Hamiltonian H = ^Z n X n Y n . The 
correlator defined in Eq. Q helps to simplify Eq. <§• 
After some easy steps, the generating function of the in¬ 
coherent FCS in the probe environment becomes: 

IF (0 = - E / % ( e_iw? - !) s Fn M (6) 


B. Renyi entropy flow 

The fluctuation relations are traditionally formulated 
in terms of entropy production that is computed using 
classical states m ■ When it comes to quantum, the 
Shannonon entropy is known to be non-linear in density 
matrix and its change is not necessarily related to the 
expectation value of any operator. This problem raises 
a careful consideration of entropy, specially that with 
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current technology developments entropy production in 
small scale systems is revealing the rich physics yet to be 
fully probed [22 , 23]. 


A generalization of Shannon entropy is the Renyi en¬ 
tropies. To evaluate the flow of Renyi entropy (R-flow) 
we need to use the perturbation theory for the M-th 
power of its density matrix. m To this end, we use a 
multi-contour Keldysh technique. We consider M copies 
of an isolated world. The contour for the degrees of free¬ 
dom of X encompasses all of the worlds and closes. This 
imposes the trace over the matrix multiplication of p x . 
For other degrees of freedom in Y, the bra and ket parts 
of the contours are closed within each world providing the 
partial trace over these degrees of freedom: that yields 
p x = Tr Y p for each world. The relevant diagrams are 
pairwise-grouped. 


The average flow of Renyi entropy during a period r 
is simply determined from Tm — (1/t) Xudt. In the 
second order we expect two interactions of the form indi¬ 
cated above Eq. ©• The two interactions can be either 
in the same world, or in different ones. The same-world 
diagrams have been considered in [TO] . The different- 
world diagrams contain contributions from quantum co¬ 
herence terms and are present when driving force is ap¬ 
plied. We studied the contribution of the quantum co¬ 
herence on the Renyi entropy flows in [24] . 



FIG. 2: (Color online) Multicontour diagrams for the R-flow 
of exponent M with interaction legs at nodes mn and n both 
in a world. Internal contours (in gray) denote system Y and 
the outer contours system X (in thermal equilibrium). Inter¬ 
actions take places at vertices m and n at different times one 
of which is fixed at t. At time t the contour of X closes in 
dotted lines. 


Single-world R-flow : The Renyi entropy flows in one- 


world diagrams are shown in Fig. §» <1) : 


M f T 

= / dt dt 'Y T *x 

JO J -oo mn 

f x ra (O (t) p %-- 1 - (t) (f) pj_ 

i ^ t p n 

x (Y m {t)Y n {t') p Y ) 

X n c t ) p x x m (t) p ™- 1 - p™x m (O x n (t) 

bn 

x (Y m {t')Y n {t)p Y ) J(ll) 


The generalised correlators in system X is defined as 


nN,M 
^mn 


(t-f) 


{x m {t')p N x X n {t)pM- N } 

Tr x W) “ 


( 12 ) 


We generalized the Kubo-Martin-Schwinger (KMS) re¬ 
lation [25] to M -worlds in [24]. The Fourier transforming 
of generalized correlator thermal in a thermal equilib¬ 
rium with temperature-independent dynamical suscepti¬ 
bility XmnM can be determined from the relation (see 
Appendix [B| : 

= exp(/3ATw)n(Mw/T)xm„(w) (13) 


Note that since the time-difference t — t' in the correla¬ 
tor is large we can shift the lower bound of time integral 
over t' from 0 — oo. 

Note that in a system with temperature-dependent 
Xmn{u) requires rescaling its temperature to T/M. This 
correlator can be easily shown to be related to the gener¬ 
alized correlators of single-world interactions S 0,M m the 
following form: S N,M (u) = exp(/3Adj)S' 0,M (cj), where 

^rnn ( u ) = Srnn (w) which is the standard spectral 
density in an environment of rescaled temperature from 
T^T* = 1 /k B fl* = l/Mfl. 

Using these definitions Eq. (11) is simplified, 


-i M r T r 1 1 , 

| lw = -— / dt dt 

‘ do J — oo 

£ / 1 (e”"- 1 ’" - 0 •«> M 

mn ' 

x (Y m (f) Y n (t) e-W-O + Y m ( t ) Y n it’) 

(14) 


Multi-world R-flow : Similarly one can calculate the dy¬ 
namics associated to multiple-world interactions. In this 
case the energy is exchanged between different worlds. 
Summing over all possible diagrams and using the gen¬ 
eralized KMS relation to simplify the result the flow of 
Renyi entropy from different-world interactions can be 
found: 
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F M \ =-~ 

I mw j- 


(Y m (: t)py) (Y n (t')p Y ) (eP* - l) 


(15) 


Details of this calculation can be found in Appendix 0 
The total flow of Renyi entropy can be obtained by 
summing over Eq. (14) and (15). By factorizing terms 
evolving with the same frequency (i.e. e ±^(4-0) the 

final result is 


t M = -m£ i)«m 

m,n J 

X (d'ran (^) H) ( 16 ) 

with A’ mn and Y mn defined in Eqs. 0 and 0, respec¬ 
tively. 

Correspondence : Comparing eq. ( p~6] ) with 0 and 0 
one can conclude the exact correspondence mentioned in 
Eq. 0. 

The case of Shannon entropy flow : Let us discuss 
here how the correspondence look like for the Shan¬ 
non entropy flow. The Shannon entropy S = Trplnp 
can be genuinely defined from the Renyi entropy in 
the following form: S' = — liniM^i SSm / dM. Using 
the correspondence 0 the flow of Shannon entropy is 
T s = (i/3) lim^o d (/$ — f c ) /df. The flow of Shannon 
entropy exactly corresponds to 


frW - 


Q 


c/n 


Q 


on 


(17) 


( 0 *) 

with Q\/ c ; the incoherent and coherent dissipated energy 
in a system of temperature T* = l/ksfl* and /3* = Mf3. 

m 


III. EXAMPLE 1: THE SIMPLEST QUANTUM 
HEAT ENGINE 

A quantum heat engine (QHE) is a system of several 
discrete quantum states connected to several environ¬ 
ments at different temperatures. The motivation for re¬ 
search in QHE comes from studying models of photocells 
and photosynthesis [26 . It has been demonstrated that 
quantum effects can dramatically change the thermody¬ 
namics of QHEs m and their fluctuations [28] manifest¬ 
ing the role of quantum coherence. 

The simplest QHE of our interest is made of a probe 
environment weakly coupled to a two level system (TLS) 
whose states are |0) and |1). The TLS itself is also cou¬ 
pled to other heat baths at different temperatures as well 
as a coherent driving force with frequency Q matching the 
two level energy difference. The specifics and the simplic¬ 
ity of the situation is that all energy exchanges take place 


by quanta Ml. The interaction between the two level sys¬ 
tem and the probe is governed by the interaction Hamil- 
tonian: H int = X 0 i (t) |0)(l|e _ * m + Xi 0 (i)|l)(0|e* nt . 



FIG. 3: (Color online) Schematics of a QHE. A quantum sys¬ 
tem with two sets of states separated by energy E\ — Eq is 
driven by external field at matching frequency. The system 
interacts with a number of environments that induce transi¬ 
tions between the states. We study the R-flows to a weakly 
coupled probe environment. 


In Ref. [24 we explicitly derived the Renyi entropy 
flow for the probe environment of this system using per¬ 
turbative expansion of the probe dynamics. Here we 
make an attempt to determine the R-flow using the full 
counting statistics method and the correspondence of Eq. 

0 . 

In the TLS, transition from an upper level to a lower 
one takes place by the operator Yio(t) = |l)(0|e zm and 
the opposite one by Yoi(t) = \0)(l\e~ int . Moreover: 

(Yio(t)Y 01 (t')p s ) = pue^-O, and <r 0 i(t)Uo(*>.> = 

Pme~ in{ - t ~ t ' ) ■ Also (Y 10 p s )(t) = p 0 ie M and (Y 01 p s )(t) = 
plop- 1 ' 11 '. Let us denote the excited and ground state 
probabilities with pi = pn and po = poo- 

Using the Kubo-Martin-Schwinger (KMS) relation [25] 
in Eq. © we can introduce the excitation transition 
rate T^ = n(Ct/T)x oi,un and the emission rate T^ = 
with the Bose function n(Q/T) = l/(e^ n — 1). 
Using Eq. 0 the full counting statistics of heat dissi¬ 
pation in the incoherent energy transfer is: 

r> (O = (*-«•" - 

( 18 ) 

Similarly, using Eq.0, the full counting statistics of 
energy transfer through quantum coherence flow becomes 


/<"•»«•) = (e-«' n -1) "(Dy (U - r t) A>i/>io 

where we used S mn , m,(- u) = e 0ul S pq , mn (u). 

Notice that Eqs. (18) and (19) are the FCS associated 
to Poisson probabilities. The reason is that the probe 
environment is weakly coupled to the quantum heat en¬ 
gine. Since we consider the dynamics to be Markovian, 
the time lag between two successive emissions in equi¬ 
librium environments at fixed temperatures is long. The 


(19) 
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events of transmissions of energy are uncorrelated. More¬ 
over due to the weak coupling the energy transfers take 
place at low transmission probability. Such a process 
can be described by Poisson probability pk = e~ n n k /k\ 
for exchanging k quanta of energy HQ, where n is the 
average number of quanta transmitted during time [0, r\. 
In the case the coupling of interaction between the two 
systems is not weak enough, or the emissions take place 
in short time intervals such that they become correlated, 
the Poissonian probability for emissions and absorptions 
are no longer valid. 

From the correspondence Eq. 0 the R-flow can be 
obtained from subtracting the two FCSs at and /?*: 

-p _ Mn(MQ/T) 

M ~ fJjM - 1 )fi/T) n(n/T) X 

— PoT^ + (r ; — r t ) poiPio) (20) 

In the second line of Eq. ( |20| ) the first two terms are 
the dissipation of heat and the third term is the energy 
transfer due to quantum coherence flow. In conclusion, 
what we calculated above matches the R-flow result we 
obtained earlier in Ref. [24] . 


IV. EXAMPLE 2: A DRIVEN HARMONIC 
OSCILLATOR COUPLED TO HEAT BATHS 


Let us consider a single harmonic oscillator of fre¬ 
quency cjo with Hamiltonian H = ujq(o) a + 1/2) is cou¬ 
pled to a number of environments at different tempera¬ 
tures with different coupling strength. We concentrate 
on a probe environment which is weakly coupled to the 
oscillator. In addition the oscillator is driven by external 
force at frequency Q. 



FIG. 4: Schematics of a harmonic oscillator of frequency uo 
interacting with environment at different temperatures. The 
environments induce transitions between the states in the har¬ 
monic oscillator. The oscillator is driven by external force at 
frequency Q. We study the R-flows to a weakly coupled probe 
environment. 


X(t)a){t) + h.c. with X being the probe reservoir op¬ 
erator. The Fourier transform of the correlator is: 
Smn(u) = f exp(—icut)Smn(t)dcu/27r. Due to conserva¬ 
tion of energy the energy exchange occurs either with 
quantum HQ or with quantum Hwq. 

We note that the time dependence of the aver¬ 
age of two operators can be written as (a)(t)d(t')) = 
((a) a))e 2u;o ^ _t ^ + (a(t)) (a^ (t')), where the (a(t)) is due 
to the driving force and therefore oscillates at frequency 
Q: (a(t )) = (a) + exp (iQt) + (a)_ exp(-iQt). This corre¬ 
sponds to the fact that the oscillator can oscillates both 
at its own frequency and at the frequency of external 
force. 

Obtaining the FCS of energy transfers is straightfor¬ 
ward from the diagrams of Fig. 0 . The incoherent and 
coherent flows are: 

+#(- Wo )((fl t a))(e“»«-l) 
+S^(ft)(a)_(a t ) + (e- in «- 1 ) 
+5(«(-0)(a) + (a t )_(e in «-l) 
-/y } (0 = 5W(fi)(a)_(at) + (e-^-l) 

Substituting these FCSs in the correspondence of 
Eq 0 using the values of and /?*, the flow of 
Renyi entropy after using the relation using S^(—uS) = 
exp ([3uj)SP (<jj) is determined to: 


= M(e^ M-1)w ° - 1) 5 (M/3 ) (w 0 ) x 

{((ato))e^-«aat))} (21) 


Given T' to be the effective temperature of the har¬ 
monic oscillator ((aaJ)) = n(uJo/T') + 1 and ((aJa)) = 
n(<uo/T f ). The KMS relation of Eq. (13) helps 
to describe the correlator in the thermal bath in 
terms of its dynamical susceptibility, i.e. = 

n(Mcj/T)x^ M ^^(cj). These help to simplify Eq. (21) into: 


frW - 

M 


Mn (Mcjq/T) x 
n((M-lH/T) u(uj 0 /T) 


{n(u 0 /r) 


n (luo/T)} 


( 22 ) 

The entropy flow is robust in the sense that it only 
depends on the probe and harmonic oscillator temper¬ 
atures and completely insensitive to external driving 
force. The entropy flow changes sign at temperature 
T = V. 


V. DISCUSSION 


We calculate the Renyi entropy flow to the probe 
environment. The coupling Hamiltonian between the 
harmonic oscillator and the probe reservoir is H(t) = 


In this paper we prove an exact correspondence be¬ 
tween the flow of Renyi ( as well as Shannon) entropy and 
the full counting statistics of energy transfers. This cor¬ 
respondence is valid for the flow to the system in thermal 
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equilibrium that is weakly coupled to an arbitrary system 
out of equilibrium subject to arbitrary time-depending 
forces. 

In the case of time-dependent external forces we need 
to introduce an auxiliary full counting statistics of en¬ 
ergy transfers. This is FCS for the case when the quan¬ 
tum forces acting on the system in thermal equilibrium 
Y ’s are replaced by their averages. The usual FCS can 
be in principle measured directly. The same applies to 
the auxiliary FCS although the measurement protocol is 
more involved. Let us describe this protocol. 

Let us notice that the forces correspond to operators 
Y n and therefore can be in principle measured directly 
as an expectation value of this observable. The output 
of this measurement is a function ( Y n (t )) which is pe¬ 
riodic with period r. From this point one can proceed 
in two ways. First way in to build an artificial system 
that interacts with system X classically and program it 
to exert classical forces on system X with values that are 
given by the results of the first measurement. One then 
collects the statistics of energy transfers to obtain the 
auxiliary FCS. The second way is more practical. One 
notices that the response of X n on the forces is linear one 
in the limit of weak coupling, so instead of measuring the 
statistics of energy transfer one can measure the matrix 
of response functions Xmn(w). Then the auxiliary FCS 
can be evaluated with the aid of Eq. (10). 

This correspondence allows us to quantify Renyi and 
Shannon entropy flows. These quantities are not acces¬ 
sible in direct measurement as they are non-linear func¬ 
tions of density matrix. Direct measurements of density 
matrix for a probe environment requires characterization 
of reduced density matrix of an infinite system, which 
is a rather non-trivial procedure and needs the complete 
and precise reinitialization of the initial density matrix. 
However, measuring the entropy flow from the R/FCS 
correspondence requires that some generating functions 
are extracted from determining statistical cumulants of 
transferred energy in experimental data. This can be 
done equally well for imaginary and real values of the 
characteristic parameter. The measurement procedures 
may be complex, yet doable and physical. 

The correspondence can have many other advantages; 
for instance: a complete understanding of entropy flows 
may help to identify the sources of fidelity loss in quan¬ 
tum communications and methods to prevent or control 
them. 

Our derivation was restricted to the second order per¬ 
turbative dynamics. There are indications the theorem 
formulated is not valid in higher orders of perturbation 
theory. It is interesting to find a similar correspondence 
that is valid in all order of interaction coupling. 
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Appendix A: A relation 

One can easily prove that in general for any multi¬ 
argument function /(cj, • • •) the following relation holds: 

J du - l) M f[w, ■■■) = 

~ [du (e^- 1 )" - l) M f(-u, ■■■) 

(Al) 

This can be easily proven by changing variable u —uj 
and using the relation between spectral density function 
of negative and positive frequencies and simplifying using 
easy algebra. 


Appendix B: Generalized KMS 

The generalized correlator of two operators A and B 
is defined (see eq. (4)): 

Sab M = J dTe™Tr{A(p)tfB(r)p^- N }/Trp^ 

This correlator in the energy eigenbasis can be rewrit¬ 
ten in matrix form 


nN,M / \ 

nm,mn v*') 


f . e -f3NE m 

J dre lUT A nm ~^pyT x 
B .^)r e- W m M 

mn Z(t3) M ~ N Z(/3M) 


2w5 (E m - E n +u) 


A R p—/3E n M 
sA'nrn - LJ, rnn c ' (3Nu> 


Z(/3M) 


(Bl) 


Z(/3) is the partition function defined as Z(J3) = 

E ie-™. 

The standard correlator is Sab (w) = 

J dre 2u;r Tr{A( 0 )H(r)p 5 }/Trp 5 becomes equal to 


Sab M = 27t<5 (. E m - E n + w) A nm B mn e l3En /Z(f3) 


where KMS relation links this to dynamical susceptibil¬ 
ity: Sab(u) = XAB(u)n(uj/T). Substituting this in (Bl) 
a generalized KMS relation is obtained: 
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S%¥ M = n{Mw/T)ef iuN XAB(u) (B2) 







Appendix C: Multiple-world dynamics 


Typical diagrams corresponding to the multi-world 
terms are listed in Fig. © 



This must be summed over all possibilities. When the 
first interaction is at the topmost world the second one 
can run between n = 2 and M. However, when we put 
the first interaction at the second topmost world the sec¬ 
ond interaction can have maximally M — 1 world distance 
with it, therefore n = 2 to M — 1. Note that we already 
consider the both positive and negative energy exchanges 
in the summation of diagrams (e-1). Extending this dis¬ 
cussion one can find the following total summation for 
all multi-world diagrams: 


M 


mw 


1 

r 




M M' 

X EE (■ Ym{t)p Y )(Y n {t')p Y } 

M'=2N=2 


X 


{ e M*-0 (5£-2 ,m H _ 2S n-i,m h + S N,M 

( S N-2,M H _ 2S n-im h + S NM j 


FIG. 5: (Color online) Typical multicontour diagrams for the 
R-flow with interaction legs at nodes m and n in two different 
worlds. 


Changing uj -> — uj in terms with indices S nrn and using 
the relation (—cj) = (cj) that can be easily 

concluded from the Fourier transforming eq. (12) and 
simplifying the summation using the KMS relation, all 
multi-world diagrams sum into 


Similar to the analysis for one world, we can write 
detailed diagrammatic values for the evolution of R-flow 
in the case two interactions occurs in two different worlds. 
After Fourier transformation the flow from the diagrams 
(e-1) of Fig. © becomes: 


M 


x (frn (t) Py ) (y n (t 1 2 3 4 5 6 7 ) Py) X 
{ e *w(t-t') [S N- 2 ,M (w) _ 2S N-1,M (w) + S N,M (w) ] 

+ e-Mt-t') [ s n-2,M (w) _ 2S n-i,m (w) + gNM (w)] j 




M 

T 


J u J -oo j rnn 

(f m (t) p Y ) (Y n (t’) p Y ) (e^ - l) 

(Cl) 


Simp lifyi ng the integration using a relation that comes 
in Eq. © we can further simplify this relation into 


^| rnw = M / (e^ M - 1)w - l) S 0 ^{u) Y mn (r,a;) 

rnn 

(C2) 
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